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^ ■ Abstract: 

H ; 

5^ , Final state qq interactions may give origin to non zero values of the off-diagonal 

element pi,-i of the helicity density matrix of vector mesons produced in e'^e~ 
annihilations, as confirmed by recent OPAL data on (p and D*'s. Predictions are 
given for of several mesons produced at large z and small pt, i-e. coUinear 
with the parent jet; the values obtained for (j) and D* are in agreement with data. 



1 - Introduction 



The spin properties of hadrons inclusively produced in high energy interactions 
are related to the fundamental properties of quarks and gluons and to their elemen- 
tary interactions in a much more subtle way than unpolarized quantities; the usual 
hadronization models - successful in predicting unpolarized cross-sections - may not 
be adequate to describe spin effects, say the fragmentation of a polarized quark. 

In Refs. m and it was pointed out how the final state interactions between 
the q and q produced in e^e~ annihilations - usually neglected, but indeed necessary 
- might give origin to non zero spin observables which would otherwise be forced 
to vanish. The off-diagonal matrix element pi-i of vector mesons may be sizeably 
different from zero [jl| due to a coherent fragmentation process which takes into 
account qq interactions; the incoherent fragmentation of a single independent quark 
leads to zero values for such off-diagonal elements. The same situation is not true 
for spin 1/2 baryons, for which the coherent fragmentation process only induces 
corrections which vanish in the limit of small transverse momentum, p^, of the quark 
inside the jet 0. Both predictions, a non zero value of pi-i for D* and possibly 
particles |^, and a value p+_ ~ {px/z^) for A {i.e., its transverse polarization) [Q 
have recently been confirmed experimentally. 

We consider here in greater details the coherent fragmentation process of qq 
produced at LEP, where the quarks are strongly polarized; we are actually able to 
give predictions for pi,_i of several vector mesons V provided they are produced in 
two jet events, carry a large momentum or energy fraction z = 2Ev/y/s, and have 
a small transverse momentum px inside the jet. Our estimates are in agreement 
with the existing data and are crucially related both to the presence of final state 
interactions and to the Standard Model couplings of the elementary e~e~^ — > qq 
interaction. 

In the next Section we review the formalism to compute the helicity density 
matrix of a hadron produced in e~ ^ qq h + X processes and give analytical 
expressions for the non diagonal matrix element pi,_i in case of final spin 1 hadrons; 
in Section 3 we obtain numerical estimates and in Section 4 we make some further 
comments and conclusions. 

2 - -i(^) the process e~e^ QQ — > V + X 

The helicity density matrix of a hadron h inclusively produced in the two jet 
event e'e^ ^ qq h + X can be written as [|l], 

i3,X,A^,A^,A^,A^,A^ 

where ^ .y y (qq) is the helicity density matrix of the qq state created in the 
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annihilation of the unpolarized e"*" and e 



Pa„,A„-;A' A'- ill) 
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(2) 



The M's are the hehcity amphtudes for the e^e^ qq process and the D's are the 
fragmentation amphtudes, i.e. the hehcity amphtudes for the process qq ^ h + X; 
the J2x,Xx stands for the phase space integration and the sum over spins of aU the 
unobserved particles, grouped into a state X. The normalization factors Nh and 
Ngg are given by 

E ^A,A,;A„A, Pa„A,-;A;,A^- (?^) ^I,A,;A;,A;- = E ' (3) 

where is the usual fragmentation function of quark q into hadron h [see also 
comment after Eq. (|22D], and 
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The center of mass hehcity amplitudes for the e qq process can be com- 

puted in the Standard Model and are given by 



^A,A,;A_A+(S, 
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+ 2g[gl (A_ + A, cos 9) - g\ g\ (cos ^ + 4A_ A 



where \fs is the total e^e~ cm. energy, Q the q production angle [i.e. the angle 
between the incoming e~ and the outgoing q) and is the quark charge. Lepton 
and quark masses have been neglected with respect to their energies and we report 
here for convenience the Standard Model coupling constants: 
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From Eqs. (0), and (|^) one finds the explicit expressions of tlie only non zero 
elements of p{qq)- 



1 - p-+;-+{.m) 
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where +, — stand for helicity +1/2 and —1/2 and where, for an arbitrary total 
energy ^/i, 
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which at = read 



a'g = 


e^ + C i9v-9A?q i9l+9l)i 


b'q = 


4C^ i9A9v)i i9v - 9a)1 


dq = 


el + e{9l-9l)q i9l+9l)i 


bq = 


^eqC9l9l 


f^q = 


2[el + e {9l+9l)i i9l+9l)q 


Vq = 


16C^ i9A9v)i i9A9v)q 


C = 




4 sin^ 9^ cos2 9^ 



Eqs. d^) and (H) hold for the production of a quark with flavour g at a cm. 
angle 9, defined as the angle between the incoming negative lepton and the outgoing 
quark; in the pt ^ limit this is the same angle as the production angle of the 
observed hadron h. However, h can be produced also in the fragmentation of an 
antiquark q and the Y.q i^i Eqs. (|I|) and takes into account also this possibility 
(g = u,d, s,c,b,u,d,s,c,b): the hehcity density matrix p{qq) for the production of 
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an antiquark at the angle 9 can be obtained from p{qq) with the simple replacements: 



P+-,-+{m,0) = pl_.^_4qq,7r-e). (11) 

The expressions @, (|^) and (p!0| ) are exact and contain both electromagnetic 
and weak interaction contributions. However, at LEP energy -y/i = M^, the weak 
contribution dominates, C ^ 1 in Eqs. (0); if one also takes into account that rjq 
is depressed by the small value of g''^ a simple approximate and useful formula for 
p_^__._^ is given by [for an exact value at -y/i = see Eqs. (^Ol)] 



Eq. ([T^p clearly shows the 6 dependence of This approximate expression is 

the same both for p{qq) and p{qq). In the case of pure electromagnetic interactions 
(i/i ^ M^) one has exactly: 

Notice that Eqs. ([l2|) and have the same angular dependence, but a different 
sign for the coefficient in front, which is negative for the Z contribution [see Eqs. 

i)]- 

By using the above equations for p{qq) [and p{qq)] into Eq. (|I|) one obtains the 
most general expression of p{h) in terms of the qq spin state and the unknown frag- 
mentation amplitudes [Q. Such expression can be greatly simplified if one considers 
the production of hadrons almost coUinear with the parent jet: the qq — >■ h + X frag- 
mentation is then essentially a cm. forward process and the unknown D amplitudes 
must satisfy the angular momentum conservation relation [Q 



where 6h is the angle between the hadron momentum, h = zq +Px, and the quark 
momentum q, that is 

sin^,~^. (15) 

The bilinear combinations of fragmentation amplitudes contributing to p{h) are then 
not suppressed by powers of {pT/{zy/s)) only if the exponents in Eq. (|T^) are zero; 
which yields 

Ax = A, - (A, - A,) = - (A; - A',) . (16) 
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In the pt ^ limit one has then the simple result for the non diagonal density 
matrix elements of spin 1 mesons |jl|, : 

Re[pi,_i(l^)] = -L^/)^„^^_z^;^ Re[p^_^_^(gg)] (17) 

x,q 

Im[pi,_i(V^)] = -^^Dio;+-^-io;-+Imb+-;-+(gg)] (18) 



with 



Nh = Y.D^,= E [l^A,A,; + j'p+-; + 4g?1 + I^A,A,;-+PP- + ;-+(?^)]- (19) 

Eq. ([T7|) and (18) explicitely show that the coherent quark fragmentation allows 
non zero off-diagonal helicity density matrix elements which, for vector mesons, 
survive also in the small pt limit; the other off-diagonal matrix elements for spin 1 
particles and all off-diagonal matrix elements for spin 1/2 particles are bound, via 
Eq. ([T^), to vanish at small pr/^/s values []1|, Recent experimental data have 
confirmed both the non zero value of pi -i{D*) |^ and the small value of p+_(A) 
§• 

In the next Section we give numerical estimates of pi,-i for several vector mesons, 
exploiting Eq. ( [T7| ) and the fact that, at least for valence quark contributions, the 
dependence on the fragmentation amphtudes either cancels out or can be expressed 
in terms of other measured quantities. 

3 - Numerical estimates of _i(^) cit ^/s = 



Let us consider Eqs. (P^ -(|19|). Despite our ignorance of the fragmentation 
amplitudes we see that in the px ^ limit, due to Eqs. (p^ ) and (|16|), only few of 
them give a leading contribution; moreover, the fragmentation is a parity conserving 
forward process, so that the fragmentation amplitudes must satisfy the relationship 



^-A,-A,;-+ = (-1)^'^+^-+^'^""- ^A,A,;+- ' (20) 

where Sf^ and Sx are respectively the spin of hadron h and of the unobserved system 
X (the intrinsic parities of the initial and final states must be the same) . In particular 
Eq. ( pO|) for spin 1 hadrons yields 

^-10;-+ = (-l)""" ^10;+- • (21) 

Notice that the parity relationship (|20|) and Eq. (0) allow to write: 

E I^A,A,;+-^ (22) 

which is the fragmentation function of quark q into hadron h, whose spin is not 
observed; such fragmentation function is independent of the quark polarization. 
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described by p{qq). Instead, the fragmentation functions of a polarized quark q into 
a hadron h with hehcity are given by: 



D 



(qq) + 1^ 



-{qq) 



X 



(23) 



which is consistent with Y.\ -Dg'^'' = D'^ and where -D'''r'' is the fragmentation 
function of quark q with hehcity into hadron h with hehcity A^^; _x -x _a (qq) 
is the probabihty for q to have hehcity A^. 

Taking into account Eq. ([T6|) and (^Ol) the above fragmentation functions read: 



D: 



E[l^lO;- 



hA 



D 



E[l^lO;- 
E[l^l2; 



0-1;+- 
h-l 



+ 1^- 



1-2; + - 



(qq) + d'^'- p- 



12;-+l 



-(99) 



h,0 



-+I 



P+-;+-(??) + 1^10;+-! P- 

+ Dg'X p-+--+{qq) . 
We now assume that, at least for valence quarks: 



(qq) 



X 

T-,/1,1 



V T^h,l 



(24) 

(25) 
(26) 

(27) 

(28) 
(29) 



The first of these assumptions simply means that quarks with hehcity 1/2 (—1/2) 
cannot fragment into vector mesons with hehcity —1 (+1). This is true for valence 
quarks assuming vector meson wave functions with no orbital angular momentum, 
like in SU{'o). The second assumption is also true in SU{Q) with = 1/2 for any 
valence q and V . Rather than taking = 1/2 we prefer to relate the value of 
to the value of Poo(^) which can be or has been measured. In fact, always in the 
^ limit, one 'has, from Eqs. (P, ([ID, (|0D, (HD and (|2|): 



Po,o{y) 



If Oq is the same for all valence quarks in V (a^ 
quark contribution: 



V 



a 



^-PoAV) 



(30) 

a^) one has, for the valence 

(31) 
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Notice that the SU (6) value = 1/2 correspond to Pq q = 1/3, that is no ahgnment, 
A = (l/2)(3poo ~ 1) = 0, for the vector meson. 

If we now use Eqs. (|l|), (||), (|2|), (|2|) into Eqs. and (H) we obtain 





1^10;+- 




- + (99) 


E,,x (1 + <) 1 


^10;+- 1 


2 



(^) = \. :^Tn '7 ■ (32) 



The numerator in the above equation depends on the squared amphtude |-Diq._,__ p 
for the — > V + X forward fragmentation process and on S^. The qq state is such 
that J = Jz = 1', the final undetected system X must then have Xx = with 
S"^ = 0, 1 or 2, the only states which can combine with the Sh = \h = ^ vector 
meson state to give a VX spin state with J = = 1. On a simple statistical basis 
these 3 possible states have respectively relative probabilities 1, 1/6 and 1/30. One 
can then conclude that the = state dominates and approximate the above 
equation ( ^2]) with 

Eg (1 + <) Dg;+ 

The actual value (^) should only be slightly smaller, due to some contribution from 



Again, if only one flavour contributes or if we can assume that does not 
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depend on the valence quark flavour, Eq. (^) further simplifies Eq. ( |5B| ) to 

ft,_.(V) |1 - p„,(l/)l ^'"^'rv^^'"^ ■ (34) 

We shall now consider some specific cases in which we expect Eq. (0) to hold; let 
us remind once more that our conclusions apply to spin 1 vector mesons produced 
in e~e~^ ^ qq ^ V + X processes in the limit of small pt and large z, i.e., to vector 
mesons produced in two jet events (e~e^ — > qq) and coUinear with one of them 
{Pt = 0), which is the jet generated by a quark which is a valence quark for the 
observed vector meson (large z). These conditions should be met more easily in the 
production of heavy vector mesons. 

Let us then start from the cases V = B*^'^, Z)*^'°. In such a case one can safely 
assume that the fragmenting quark is the heavy one so that Eq. (^) applies and 
one has: 

Pi,_i(i?*+) [l-Po,o(5*+)]P+-;-+W 

p,,_AB*-) ^ [l-Po,o(5*")]P+-;-+W (35) 
Pi,-i(5*°) ^ [l-Po,o(5*°)]P+-;-+(^&) 



Pi,_i(D*+) :^ [l-po,o(^*^)]P+-;-+(cc) 

p,,_,{D*-) ^ [l-po,o(^*")]P+-;-+(cc) (36) 

Pi,-i(^*°) ^ [l-Po,o(^*°)]P+-;-+(cc) 
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Similarly one obtains: 

Pi,-i(0) - ^ [1 - Pofii^)] [P+-,-+iss) + p+_;_+(ss)] (37) 

where we have assumed Df^l_ = Dp^]_, as it should be. 

For p's, assuming all valence quark fragmentation functions to be the same, one 

has 

Pi,-i(p^) - ^ [1 - Po,o(P^)] [p+-,-+iuu) + p+-,-+{dd)] 

Pi,-i(P°) - ^ [1 - Po,o(P°)] [p+-,-+iuu) + p+-.-+{dd) 

+ P+--+iuu) + p+--+idd)] (38) 
Pi,-i(p") - ^ [1 - Po,o(P")] [p+-,-+idd) + P+-,-+{uu)] . 

The assumption that all valence quark fragmentation functions are the same is 
very natural for p's, but it might be weaker for K* mesons; if nevertheless we assume 
that, at least at large z, Df_^_ = -D^+ and similarly for K*^ and K*~ , we have 

Pl,-liK*'-) ^ i[l-po,o(^*+)][P+-;-+W+P+-;-+(^"^)] 

Pi,-i(^*°) ^ ^ [1 - Po,o(^*°)] [P+-,-Add) + P+-;-+(s-^)] (39) 

Pi,-iiKn ^ ^[l-p,^,iK*-)][p+-,-4^u)+p^-,-4ss)]. 

A predominant contribution of the s quark would instead lead to results similar to 
those found for B*. 

Eqs. (|35D-(p9D show how the value of Pi^i{V) are simply related to the off- 
diagonal matrix element p+---+{<l<i) of the qq pair created in the elementary e~e^ — >■ 
qq process; such off-diagonal elements would not appear in the incoherent indepen- 
dent fragmentation of a single quark, yielding p^ _i(^) = 0. 

We can now make numerical predictions by inserting into the above equations 
the exphcit values of p+- -+{qq) at ^/s = M^, Eqs. (||), ([ToD and (^ with sin^6';^ = 
0.2237, = 91.19 GeV^ = 2.50 GeV §: 

sin^^ 



P+-; 


-+iuu) = 


-0.361 


[1 -0.013i 


P+-; 


-+{uu) = 


-0.361 


[1 + 0.013i 


P+-: 


■_+{dd) = 


-0.171 


[1-O.OlOi 


P+-: 


■^^{dd) = 


-0.171 


[1 + O.OlOi 



;i + cos2^) 0.29 cos ^ 

;i + cos2^) -0.29 cos 

1 + cos2e) + 0.39 cos e 
sin^ 9 

;i + cos2^) -0.39 cos ^ 



(40) 
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The values for s, b and c quarks are respectively the same as for d and u. 

If we instead use for simplicity the approximate expressions (p!2D, valid at y/s = 
and which are the same for p+_;_+(gg) and p^- -^{qq), we have the simple 
results 







, V -, Dili (7 

-0,170 |l-P.,«(B-)li^,„,., 








, -. olli (7 

- 0.109 ±0.015 — 

1 + COS"' i) 


(41) 






V / J, N T bill (7 

-0.360 ll-Po.o(/>-)li^,„,., 








-(0.216 ±0.007) ' 

1 + COS"' y 


(42) 


Pi,-i(0) 




r /NT Sin (7 

-0,170 |l-Po,oWli^,„,,, 








sin^ (9 

(0.078 ±0.014) 

^ ' l + cos2^ 


(43) 


Pi.-i(p^'°) 




sin^ Q 

-0,265 |l-Po,o(p)li^,„,,, 


(44) 


Pi,-i(^*^) 




sin^ ^ 

-0.265 [1 - Mi.*)] ^^^^^^^^ 


(45) 


Pi,-i(^*°) 




-0.170il-Mir*)]^-J,^ 


(46) 



where we have used Po,o(^*^'°) = 0.36±0.09, Po,o(^*^'°) = 0.40±0.02 and Po,o(</') = 
0.54 ± 0.08 no data are available on Poo(p) ^-iid pQQ(i^*). Notice that in such 
approximation p^ _i(^) is real and that the cos Q term in the denominator of Eq. (§) 
has been neglected. This term would induce small differences between the values 
of pi _i(i?*+) [or pi _i(Z}*+)] and p^_^{B*^) [or p^_-^(p*^y\\ it has much smaller 
effects on the values of _i(0), Pi -i(p) and p^ ^^(iC*). 

Finally, in case one collects all meson produced at different angles in the full 
available Q range (say a <Q < it— a, \ cos6\ < cos a) an average should be taken in 6, 
weighting the different values of p^ _i(6') with the cross-section for the e~e~^ V+X 
process; this amounts essentially to weight the values of p+---+{qq', 0) appearing in 
Eqs. (P3D-(PUD and given in Eq. ( ^0]) or ([T^ with the cross-section for the e~e"'" qq 
process, proportional to the normalization factor A^^^^ given in Eq. (H). Such an 
average has a simple analytical expression if one uses the approximate value ([I^): 

(Pi,-i(5*^'°))K.-a] ^ -(0.109 ± 0.015) (47) 

3 + cos"' a 

(Pi,-iP*^'°))k.-.] - -(0.216 ± 0.007) (48) 
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(Pi-i(0))[a,.-.] ^ -(0.078 ± 0.014) ^-^^ (49) 

o + cos^ a 

(Pi,-i(p^'°))K.-a] ^ -0.265 [l-po,o(p)] 77— (50) 
(Pi,-i(^*^))[...-«] ^ -0.265 [l-po,o(^*)] 77— (51) 

^ -0.170 [l-Po.o(^l o^,,,2^ • (52) 

We have explicitely checked that the full expression (^OD yields almost identical 
results [and a negligible imaginary part]. 

4 - Comments and conclusions 

We have computed, within a general factorization scheme, the off-diagonal helic- 
ity density matrix element pi,-i of vector mesons produced in e~e+ —>■ qq V + X 
annihilation processes; such element can be - and in few cases has been - measured 
via the angular distribution of two body decays of the meson in its helicity rest 
frame. Our results hold for small pt and large z hadrons, in particular we expect 
them to hold for heavy mesons which should more easily satisfy such requirements. 

Our results for 0, Eq. (^), are in agreement with data, Repj^__]^(0) = — 0.11±0.07 
P]; notice that such data refer to values of 2; > 0.7 and cos a = 0.9, but still have 
large errors. Our results for D*, Eq. (|48|), have the same negative sign, but are 
larger in magnitude than the value found by the OPAL collaboration, Rep^^ _i{D*) = 
—0.039 ± 0.016 0. There are good reasons for that: data on D* are collected for 
z > 0.5, and might still contain events to which our calculations do not apply and 
for which one expects p^ _i = 0; one should also not forget that our predictions are 
somewhat lessened (in magnitude) by contributions from 5*^ = 1 [see comment after 
Eq. (H)]. 

We notice that while the mere fact that pi,_i differs from zero is due to a coherent 
fragmentation of the qq pair, the actual numerical values depend on the Standard 
Model coupling constants; for example, pi,-i would be positive at smaller energies, 
at which the one gamma exchange dominates, while it is negative at LEP energy 
where the one Z exchange dominates. pi,-i has also a peculiar dependence on the 
meson production angle, being small at small and large angles and maximum at 
e = n/2. 

Such coherent effects in the fragmentation of quarks might not play a role in 
unpolarized observables, where they are usually neglected; however, they should be 
taken into account when dealing with more subtle quantities like off-diagonal spin 
density matrix elements. Many of these effects vanish in the limit of small intrinsic 
momentum of the hadron inside the jet, px/Eh — * 0; this happens, for example, 
in the fragmentation of quarks into spin 1/2 hadrons [0]. The quantity considered 
here, instead, survives also in the small pt limit; we actually exploit such a limit in 
order to make numerical predictions. 
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The recent data are encouraging; it would be interesting to have more and 
more detailed data, possibly with a selection of final hadrons with the required 
features for our results to hold. A measurement of the pt of final hadrons and a 
study of the dependence of several observables on its value would offer many more 
possibilities of testing both the dynamics of the fragmentation process and unusual 
aspects of the basic interactions; a measurement of p+_(A) with a selection of A 
particles with 7^ is already available 0] and in agreement with our expectations. 
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